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Abstract 

A scanning pencil -beam scatterometer for ocean windvector determination has 
potential advantages over the fan-beam systems used and proposed heretofore. The 
pencil beam permits use of lower transmitter power, and at the same time allows 
concurrent use of the reflector by a radiometer to correct for atmospheric 
attenuation and other radiometers for other purposes. The use of dual beams based 
on the same scanning reflector permits four looks at each cell on the surface, 
thereby improving accuracy and allowing alias removal. 

This paper describes simulation results f or a spaceborne dual-beam scanning 
scatterometer with a 1-watt radiated power at an orbital altitude of 900 km. Two 
novel algorithms for removing the aliases in the vindvector are escribed, in 
addition to an adaptation of the conventional maximum- I ikel ihood algorithm. The 
new algorithms are more effective at alias removal than the conventional one. 
Measurement errors for the wind speed, assuming perfect alias removal, were found 
to be less than 10%. 

1.0 INTRODUCTION 

Plans for wind-vector measurement with spaceborne radar scatterometers often 
call for measurements with a radiometer on the same spacecraft. Although many of 
these measurements may be independent of the scatterometer , at least some of them 
should be used to correct errors in the measured barkscatter due to atmospheric 
attenuation. In fact, a radiometer should always be used in conjunction with a 
scatterometer for this purpose. 

A major problem with this correction occurs when the radiometer and 
scatterometer have different scan patterns and the radiometer has a larger 
footprint than the scatterometer [Moore, et al . , 1983], The problem exists because 
storm cells are often small compared with the size of a single radiometer 
footprint, so combining several radiometer measurements to correct a scatterometer 
measurement for a cell overlapping several radiometer cells results in significant 
errors in the "correction". The instrument configuration discussed here is 
intended to overcome this problem. 

This arrangement calls for both instruments to have the same scan pattern and 
for the radiometer to have a footprint coincident with that of the scatterometer 
and, to the extent possible, a footprint of the same or smaller size. This is 
achieved by using the same scanning pencil -beam antenna for both scatterometer and 
radiometer. In addition, a second beam achieved by using an offset feed in the 
same antenna allows the scatterometer to have four rather than two azimuth angles 
(relative to the wind direction) for viewing each scatterometer cell. 

The basic configuration is shown in Figure 1. Two circles at different 
distances from the suborbital track are shown, along with a line parallel to the 
track but displaced from it. The outer circle is the locus of scan positions for 
the radiometer and for scatterometer beam no.l, and the inner circle is the locus 
for scatterometer beam no. 2. The intersections of these circles with the 
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Fig. 1: Basic Configuration of 

Scanning Scatterometer 


Fig. 2: Flow Chart for 

Simulation 


displaced track line show the four different angles that are used by the 
scatterometer to view each cell on the surface. Use of only two look directions 
results in up to four "aliases" of the wind direction [Wurtele, et al., 1932] Use 
of additional beams can reduce this alias problem and allow determination of the 
correct wind direction much of the time without use of collateral data or pattern 
recognition [Shanmugan, et al., 1982] This paper presents simulations that 
illustrate how well the four-beam system solves the alias problem. 

Advantages of this system include: 

(1) Only one antenna is needed for both radiometer and scatterometer. 

(2) Use of the relatively large aperture required for the radiometer allows 
the scatterometer power to be quite low (1 watt in these simulations). 

(3) Since all scatterometer measurements are at only two fixed angles of 
incidence, errors in modeling the incidence-angle variation of the 
scattering coefficient do not exist, and the data processing for the 
scatterometer need not have algorithms to handle incidence-angle 
variations. 

The primary disadvantage of this system is that the look directions for the 
different cells are neither orthogonal (as on SEASAT) nor constant. On the other 
hand, the effect of the rotation of the earth on the SEASAT observations was to 
make the orthogonal beams on the spacecraft more difficult to handle than they 
would have been if they had been oriented relative to the ground track rather than 
relative to the orbit plane. Hence the complexity in the computations associated 
with non-orthogonal beams is at least partly cancelled by the complexity 
associated with correcting for earth-rotation effects in the orthogonal -beam 
system. For some angular combinations with the scanning beam the accuracy is less 
than with orthogonal looks, but these may be avoided at some sacrifice in swath 
width. 

Here we treat two subjects associated with the scanning-beam system: 
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algorithm development for determining wind vectors with the scanning four-beam 
scatterometer and measurement errors with such a system. The treatment is 
accomplished using Monte Carlo simulations. All simulations assume an effective 
transmitter average power of 1 watt at a spacecraft altitude of 900 km. Other 
parameters used are listed in Table 1. The algorithms used include a modification 
of the maximum-likelihood (or SOS) algorithm used for SEASAT [Jones, et a 1 . , !d.° 
and two simple pattern-recogni tion schemes based on comparison of results from 
nearby cells. The averages of norms between the true wind vectors and the 
simulated ones are computed to illustrate the errors. 

2.0 SIMULATION METHOD 


The major steps in the simulation of the scanning scatteroineters are shown m 
Figure 2. 

Step 1: Input o f System Paramet ers: This step inputs all parameters that 

are needed in the following computation. The list of the parameters is shown in 
Table 1. 


Step 2: Compute Noise-Free Backscatteri nq Coefficients {o*J} : This step 

computes a set of noise-free backscatteri ng coefficients in equation (1). 




(d8) = 10 [G(* q - 


180 °.°r c t ) + 


100°, 0 , e ) log 


10 


iu 0 |] 


( 1 ) 


where : 

0 = true wind vector 

s W 0 \cos* 0 ,\Q 0 \*\n* 0 ) 

; = true wind direction relative to suborbital track (see Figure 1) 

■'ri = look angle relative to suborbital track (see Figure 1) 
i = notation to distinguish looks 
- incidence angle 
ci = polarization (HH or VV) 

G(-),H(-) = G-H Table [Schroeder, et al., 1982] describing the scattering 
model function 


TABLE 1 

PARAMETERS OF SIMULATION 


Transmit Power (W) 

1.0 

Noise Temperature (°K) 

273.0 

System Loss (dB) 

4.0 

Noise Figure (dC) 

4.0 

Seamwidth-H (deg) 

1.0 

Beanwidth-Y (deg) 

1.0 

Radiation Efficiency 

0.8 

Fraction of time to pass footprint kj 

0.1 

Signal and Noise Integration Time Ratio k^ 

10.0 

Wavelength (MM) ' 

22.0 

Satellite Height (km) 

900.0 

Slant Range (km) 

1127.6 


TO 
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Step 3: _ Compute flornial i z_e(l_St^nddrdJ)eviati on of Measurement rrror Kp(t) Pin? 
t o Poi se. This *> L < * t > computes the normalized standard deviation Kp(t) of 
measurement error [Grantham, ot al.,] as 


K C (i.) = 


■1 2 -2 
L + SNR -_ 
j , r T irr . k ., 

cl c I N 


>NR 


( 2 ) 


where: 

SNa = signal -to-noise power ratio 
i' C = receiver bandwidth 
T[ = signal integration time 
kp = noise-to-signal integration time ratio 
i = notation to distinguish looks. 

Note that only the receiver noise is taken account of in these simulations, 
although the measurements might also be noisy because of sampling variability. 


S tep 4 : Co mpute il oisy flac ksca tt ering Coeff icients j o%). This step computes 

the set of noisy back sea 1 t er Tmj”coe f f i c i ent s { o9] in equation (3). It is assumed 
that arc distributed log-normally. 

In : In o° e + Ao° (3) 

where : 

Ao 9 = Gaussian random variable with 0 mean and K (t) deviation, 
p = notation to distinguish looks. 

Step 5: Compute Hin d Alias from |o°} . The maximum-likelihood method is 

applied as the algorithm for computation of the wind alias from {o^} . Thus, the 
wind vectors !)( = ( j Oj cos<j>, | fl| sin<t») ) at which the probability density function (PDF) 
of equation (4) have local maxima are computed [Jones, et al., 19H2] 

exp { - l (lno°-f(U,t)) 2 /2K |i 2 (4)} 

PDF ( U ! i a 3 1 } ) = — - 4 - (4) 

/ dll exp 1 - E (lna°-f(il,£)) 2 /2K 2 (t)} 
i-\ * 

whore: 

f(0,4) = (lnlO)[G(»-iP t - 18O° f 0 lf c £ ) + H(t-<l< r lfl0 0 ,O rejl )].log 10 fu| (5) 

The procedures from Step 1 to Step 5 are iterated for several true wind 
vectors ^ ground points, incidence angles and two polarizations, because these 
parameters seriously influence wind alias and wind measurement error. Other 
parameters are fixed at suitable values for spaceborne systems. Wind-alias 
removal and measurement error are estimated in Steps 6 and 7. 


3.0 WHIP- ALI A S REMO VAL 

The PDF of equation (4) usually results in two or four local maxima of the 
wind direction, resulting in wind aliases. The local probability maximum for a 
wrong wind vector often has a larger value than the local maximum of the correct 
wind vector. Hence, it is important to remove wind aliases by other means. 

ixamples of the simulation results are shown in Figures 4-13. Parameters of 
the simulation are shown in these figures and in Table 1. In these figures 
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"ground-track " shows the direction of spacecraf t motion and "A7.Ano." shows 
f orw a rd azimuth look angles ol oath boom. ( Ihe reference axis is the ground- 
track IT Aft azimuth look angles of each beam are given by subtracting forward 
look angles trail l.'‘»0°. I he angle between the vector and "grounn-track" shows the 
wind direction and the magnitude of the vector shows the wind speed. 

'lote that the type of "box" used for simulating the pattern recognition 
techniques is not completely realistic because it is based on the angles used for 
the simulations. The horizontal (j) dimension on the "chicken tracks" corresponds 
with along-track spacings for the satellite. The vertical (i) dimension on the 
"chicken tracks" corresponds to across-track spacings for the beam, but these 
dimensions are based on the angles selected for simulation rather than on uniform 
spacings. Figure 3 illustrates the dimensions of these boxes for different 
central angles chosen. A more realistic simulation would involve square boxes, 
but would mean repeating the simu 1 itinns for these angles. 


A • . • n L . k A. ij> (T* greel 
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Fig. 3: Ground Spacings of Cells vs Azimuth Angle, (a) Outer Cells 

for Window 1; (b) Outer Cells for Window 2 


The 4-vector plots ("sea chickens") shown in Figure 4a were prepared for each 
case studied, but only one example is presented here. These show the wind alias 

vectors, and X shows any wind al ias whose normalized probability P(K), is 

over 0.1, where 


p{k) - 

n 

: PDF ( K) 

K-l 

where : 

PDF = probability density function in equation (4). 

K = notation to distinguish alias. 

K n * number of aliases. 

The "Alias Removal" code shows the wind vectors which are 
following alias removal algorithms. 


( 6 ) 


selected by one of the 


( 1 ) Algorithm 1 - (Maxi mum Probability) 

tiU.J.Kj «- Max[?(I,.J,K) j 

c K=K K 
c 


( 7 ) 
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In this algorithm, the alias that has maximum normal i 2 nd probability is 
selected as the correct wind vector. Where: 

!'( I ,.) ,K 1 - wi nd al i as , 

I ,J : notation to distinguish ground point (see figures 3(a) and (b)), 
K - notation to distinguish alias, 

?(I,.J,K) = normalized probability of wind alias. 


r i -j-j re 4b shows an example of the wind vectors obtained by selecting the 
alias in Figure 4a using Algorithm 1 for 4 m/sec winds with vertical polarization 
for the forward beam pointed in the upwind direction. Rince the results with 
Algorithm 1 are generally inferior to those with the other algorithms, the results 
obtained witn this algorithm are not repeated for the other examples shown. 


( 2 ) Algori t h m 2 - (Patte rn Recog nition - 1) 


*J(X,J.K C ) 


K=K 


Min ; C J Mi n { | > J ( i , j ,k ) - l3( I,.),K)!/|IJ( I.J.K) | } ] ] 
ij k 


(P) 


In this algorithm, the wind alias v/hich is satisfied with equation (V 
selected as the correct wind vector. There T[*] is the average over i and j and 


the ranges 

of i and j 

are as follows: 




(i) 1-1 

s i s 1+1, 

j = J-l, J 

= J+l 

-Window 1 

(9) 

(ii) I-l 

< i < 1+1, 

j = J-2, j 

= J+2 

-Window 2 

(in) 


(See Figure 3) 


Window 1 uses adjacent cells in the pattern-recogr 1 tion scheme, whereas Window 2 
uses both adjacent cells and the next cell away in the cross-track direction. 

Thus, Window 1 uses 9 cells for alias removal, whereas Window 2 uses IS cells. 

This is particularly important for the large azimuth angles, where the geometry of 
the simulation indicated in Figure 2 causes adjacent cross-track cells to be very 
close together. Of course, tne assumption of uniform wind fields is more 
important in using Window 2. 

Window 1 is ised in Figures 4c and 5a to illustrate the use o. Algorithm 2 
(Pattern Recognition - 1). The other examples shown are for Algorithm 3 (Pattern 
Recognition - 2), since it gives better results. 


( 3 ) Algorithm 3 - (Pattern Recognition - 2) 


U( I ,J,K ) « Min[E [Min }-*n{P( I ,J ,K) -P(i , j ,k ) } • } 

c K*K K ij k 
c 

•!$(1,j,k) - U(I,J,K)|/|I?{I,J.K)| }]] 


(ID 


In this algorithm, the wind alias which is selected by equation (11) is 
presumed to be the correct wind vector. The range of i,j is given by equation (9) 
or (10). This weighted algorithm gives generally better results than Algorithm 
2. Therefore it is used in all of the examples following the initial ones that 
are shown in Figures 4 and 5 to illustrate the improvement of Algorithm 3 over 
Algorithm 2. Window 2 is used in Figures 6b, 7b, 10 and 13, since it gives 
significant improvement when the forward beam is in the crosswind (90°) direction. 
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The simulation results show the following: 

(1) F.very algorithm is apt to select wrong wind vectors in the case of low 
wind speed. The error pattern of algorithm 1 seems to be random, but 
the error patterns of algorithm 2 and of algorithm 3 seem t> be meru 
consistent. Thus, algorithm 2 and algorithm 3 are apt to select wind 
vectors 180° from the correct direction. This error can easily b,> 
removed by examining the wind patterns. 

(2) Algorithm 2 is apt to select opposite wind vectors even in the case of 
high wind speed. 

(3) Window 2 iinptoves the performance of Algorithm 2 and of Algorithm 3. 

(4) Algorithm 3 with window 2 will select correct wind vectors almost 
perfectly. 

Note that the low-wind-speed errors might be reduced by use of higher power than 
the 1 watt assumed. They would also be smaller if th' satellite were at a lower 
altitude than 900 km. 


4.0 MEASUREMENT ERROR Of WIND VEC TOR 


In this section, the measurement error of the wind vector (after alias 
removal) is investigated. It is presumed that the al ias-removal algorithm works 
perfectly. Scanning scatterometers with one beam (two look angles) are also 
investigated here. 

The simulation results are shown in Figures 14-16. Parameters of the 
simulation are shown in Table 1 and in these figures. In these figum the 
horizontal axes show the forward look angle of the inner beam. The vertical axes 
show measurement error MC 


MC 


where: 


^[|u| - !iJ 0 | 2 ] 


(i?) 


0 o = true wind vector. 

0 = simulated wind vector. 

Twenty-five noisy samples of wind vectors were simulated in look directions 
relative to the wind of 0°, 30°, 60° and 90°. Hence, one hundred samples are used 
in the average of equation (12). These figures show that vertical polarization, 
small incidence angle, and high wind speed ?educe measurement error of 
scatterometers because they increase received power. Higher power at larger 
incidence angles might improve performance, but no such simulations were 
performeJ. 

The averages of 2400 simultion samples as to three wind speeds (4.0 m/s, 12.0 
m/s, 24 m/s), four wind directions (0°, 30°, 60°, 90°) and eight azimuth look 
angles (10°, 20°, 30°, 40°, 50°, 60°, 70°, 80°) are shown in fable 2. The rows of 
Table 2 show parameters (Incidence angle and polarization) of the inner beam and 
the columns show parameters of the ou*:er beam. Although a smaller incidence angle 
is better for scatterometers, scatterometers whose Incidence angle is 50° for the 
outer beams were Investigated because 50° is o*ten used in radiometers, and 
because 50° gives a wider measurement swath. These figures and Table ? obviously 
show the following: 

(1) Two beam (four-look) scatterometers are better than one-beam (two-look) 
scatterometers. 

(2) Two-beam scatterometers give quite low reasurement error over look 
angles from 20 0 to 80° off the surface track. 
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(3) The simulation results in the case of low wind speed are not good, but 
an increase (i.e., to two watts) in transit, t power would resolve this 
problem easily. 

TABLE 2 

AVERAGE MEASUREMENT ERROR 
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5.0 CONCLUSIONS 

Wind alias removal and measurement error of wind vector in scannlno two-beam 
scatterometers were estimated using computer simulation. In the estimation of 
wind alias removal, two novel algorithms using pattern recognition technique • re 
introduced and simulation results showed these algorithms are more powerful tt.an 
the conventional maximum-likelihood algorithm. Note that these algorithms would 
not be powerful in scatterometers with fixed fan beams like StASAT SASS. In the 
estimation of measurement error of wind vector, simulation results showed the 
scanning scatterometers give low enough measurement error over an arc from 20° to 
the side of the spacecraft track to 80° to the side to be practical. The average 
of measurement error was under 10.0%. 

Only receiver noise was taken account of as the primary factor of noise in 
this simulation. Other noise sources should be taken account of in future 
studies. Additional noise would increase K p . However, it is possible to 
reduce Kp by increasing transmit power. Only 1-watt transmit power was used in 
this simulation. More detailed studies would make certain that the scanning 
scatterometers are very powerful spaceborne wind-vector-measurement systems. 
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Figure 4: Examples for Different Stages In Allas Removal. Wind Speed 4 m/sec. 
Look Angle 0° (Upwind for Forward Beam). Vertical polarization. 

(a) Result of basic algorithm, showing aliases; (b) Aliases with 
maximum probability on basic algorithm; (c) Aliases selected by 
algorithm 2; (d) Aliases selected fay algorithm 3. 
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Figure 5: Examples for Algorithms 2 and 3. Wind Speed 4 p/sec. Look Angle 

30°. Vertical Polarization, (a) Algorithm 2; (b) Algorithm 3. 
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Figure 6: Examples for Algorithm 3, Windows 1 and 2. Wind Speed 4 m/sec. Look 
Anglo GO*. Vortical Polarization, (a) Window Is (b) Window 2. 
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Figure 7: Examples for Algorithm 3, Windows 1 and 2. Wind Speed 4 m/sec. Look 
Angle 90°. Vertical Polarization. This look angle is always the 
worst case, (a) Window 1; (b) Window 2. 
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Figure 8: trample for Algorithm 3, 
Window 1. Wind Speed 12 
m/s. Look Angle 0*. 
too Horizontti Polarization. 


Figure 9: Example for Algorithm 3 
Window 1. Wind Speed 12 
m/s. Look Angle 60* 
Horizontal Polarization 
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Figure 10: Example for Algorithm 3, 
Window 2. Wind Speed 12 
m/s. Look Angle 90°. 
Horizontal Polarization. 


Figure 11: Example for Algorithm 3, 
Window 1. Wind Speed 24 
m/s. Look Angle 0°. 
Horizontal Polarization. 
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Figure 12: Exanple for AlgorlthR 3, 
Window 1. Wind Speed 24 
m/s. Look Angle 60*. 
Horizontal Polarization. 
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Figure 13: Example for Algorithm 3, 
Window 2. Wind Speed 24 
m/%, Look Angle 90*. 
Horizontal Polarization. 
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Figure 15: Errors for Kind Speed of 12 m/sec vs Azimuth Angle with Beam Incidence 
Angles and Polarizations as Parameters. 
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Figure 16: Errors for Mind Speed of 24 e/sec vs Azlnuth Angle with Bean Incidence 
Angles r*d °o1ar1zat1ons as Paraneters. 
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